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ABSTRACT 



Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming 
regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us 
to study the star formation efficiency at sub-galactic scales around redshift z = 1.2. We present a method to obtain the gas and star 
formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding 
scales are not resolved. This method is based on the identification of these structures in position- velocity diagrams corresponding to 
slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star 
forming distant galaxies - EGS13003805, EGS13004291, EGS12007881 and EGS13019128 in the AEGIS terminology - to determine 
the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity 
is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [Oil] recombination line. We identify 
the ensembles of clumps with the angular resolution available in both CO and [Oil] spectroscopy, i.e. 1-1.5". SFR and gas surface 
densities, E S fr and S gas , are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM 
slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation at a scale of ~ 8 kpc. The data points globally follow a 
power law E S fr £g as of exponent N=l corresponding to a depletion time of 1.0 Gyr, but with a large scatter, which means that the 
depletion time varies from point to point within the galaxies. We find a depletion time of 1.4+0.9 Gyr. 

Key words, galaxies: evolution - galaxies: high redshift - galaxies: ISM - galaxies: starburst 



1. Introduction 

Ten billion years ago, between redshift 1 and 3, observed galax- 
ies formed their stars at rates as high as ten times that o f 
the Milky Way today (INoeske et al. 20071 IDaddi et al. 20071) . 
This implies a more abundent gas supply, either fueled 
by major mergers or by a semi-continuous gas accretion 
dDekel et al. 2009l|Dekel, Sari & Ceverino 2009| ). Recent obser- 
vations of normal, massive, star-forming distant galaxies tend 
to show that most of these galaxies O 80%) are not ex - 
periencing major mergers or interactions (ITacconi et al. 20131) , 
and that their light profiles are similar to those of their low 
redshift counterparts ( |Wuyts et al. 2011b) . Also, the Kennicutt- 
Schmidt (KS) relation between total molecular gas and star for- 
mation rate_(SFR)_surfaceto 

shifts dKennicutt 1998al |Wyder et al. 20091 IGenzel et al. 20101 
ITacconi et al. 2010k implying that the star formation processes 
seem to be largely independent of the cosmic epoch. The aver- 
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age fraction of molecular gas relative to the total baryonic mass 
must therefore have been up to ten times more at high redshift 
than in today's nearby galaxies, which is indeed corroborated by 
direct observations (ITacconi et al. 20131 . 

Numerical simulations suggest that typical massive star 
forming galaxies at high redshift continually accrete gas from 
the inter-galactic medium along cold and clumpy streams stem- 
ming from the cosmic web, while the disk fragments into 
a few clumps dKeres et al. 20051 Bournaud & Elmegreen 2009 , 
IDekel et al. 20091 IVan de Voort et al. 20111) . Galactic disks of 
distant galaxies are indeed found to be fragmented in a num- 
ber of giant molecular clouds, or clumps, which differs from lo- 
cal galaxies, in which the molecular gas is scattered in numer- 
ous lower mass giant molecular clouds. The high-redshift star- 
forming complexes are found to have typical scales of ~1 kpc, 
masses u p to 10 9 M , and contribute to 10-25% of the galaxy lu- 
minosity (Forster Schreiber e t al. 20111) . Most studies only con- 
sider averaged quantities for distant galaxies, as direct observa- 
tions of the star-forming regions of high redshift galaxies are 
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Fig. 1. Composite HST images of the four galaxies studied here, 
combining ACS I and V bands. The green box shows the DEEP2 
slit used for each galaxy, whose size is indicated in Table [T] 
(Appendix). We assumed that the DEEP2 survey measured the 
position of the center of the galaxy more accurately than the 
HST, and shifted the HST image accordingly. North is up, east 
is left. 



still challenging. KS relations have been derived at sub-galactic 
scales for nearby galaxies ( [Bigiel et al. 2008| ) but not yet exten- 
sively at high redshift. 

In this letter, we use IRAM Plateau de Bure CO observations 
from the PHIBSS survey (ITacconi et al. 2013b and Keck DEEP2 
spectra (Newm an et al. 2012b of four massive galaxies at red- 
shift z ~ 1.2 in order to investigate the star formation efficiency 
within their clumps, or groups of clumps. More precisely, our 
goal is to determine the SFR and the gas mass surface densities 
at sub-galactic scales and examine their correlation. 

The four galaxies were drawn from the All-wavelength 
Extended Groth strip International Survey (AEGIS), which pro- 
vides deep imaging in all major wave bands from X-rays to 
radio, including Hubble Space Telescope (HST) images, and 
optical spectroscopy (DEEP2/Keck) over a large area of sky 
(INoeske et al. 20071 iDavis et al. 20071|Cooper et al. 20TT| 2012, 
INewman et al. 2012b , thus providing a complete set of galaxies 
between 0.2 < z < 1.2. The four galaxies studied here were 
selected from the sample analyzed by ITacconi et al. 20101 and 
2013, which gathered non major-merger, luminous, star forming 
galaxies at z~1.2, with stellar masses above 3. 10 10 M© and SFR 
above 40 M© yr" 1 . These four galaxies - named EGS 13003 805, 
EGS 13004291, EGS12007881 and EGS13019128 in the AEGIS 
terminology - correspond to the massive end of the 'normal' star 
forming galaxies at z~1.2 (ITacconi et al. 2010"b . Figure [T] shows 
them in the I and V band of the HST Advanced Camera for 
Surveys (ACS), as well as the DEEP2 slits we used. To compute 
the physical distances, we adopt a cosmology with Q A = 0.73, 
£l m = 0.27, and H = 70 km.s^.Mpc -1 . 



2. Determination of the total gas mass and the SFR 

2.1. The gas mass 

We use high resolution observations of the CO(3-2) transition 
shifted into the 2 mm band for z~ 1 .2 sources carried at the IRAM 
Plateau de Bure, including the most extended 'A configuration. 
The angular resolution obtained is between 0.5" and 1.5" for 
all galaxies (Table [T| in Appendix), corresponding to physical 
scales between 4 ana 13 kpc. The CO luminosity associated to 
any region of flux can be expressed as 



K.km.s .pc 2 



3.25 10 7 (1 +z)~ 



Sco^V 
Jy.km.s -1 



/ Yobs V 2 I Dl \ 
I GHz/ \Mpcj 



where ScoAV is the velocity integrated fl ux, v b s the observed 
frequency and D L the luminosity distance ( Solo mon et al. 1997b . 
To derive H 2 m asses, we follow the approach used by 
ITacconi et al. 20101 This assumes a Milky Way like conversion 
factor of a = 3.2 M (Kkms -1 pc 2 )" 1 between the CO(l-O) 
luminosity and the H 2 mass, a factor 1.36 to account for in- 
terstellar helium, and a further correction by a factor 2 for the 
CO(3-2)/CO(l-0) luminosity ratio. The total gas masses of the 
four studied galaxies are indicated in Table [T] in the Appendix. 

2.2. The star formation rate 

Ha or [Oil] recombination line luminosities are direct and sen- 
sitive tracers of young massive stars, and thus of the SFR 
(IKennicutt 1998bb . Since the Ha line from the four galaxies 
considered here lies in the middle of a low atmospheric trans- 
mission band, ground-based Ha spectroscopy would be impos- 
sible, so we estimate the SFR from the [Oil] recombination 
line (rest frame wavelength 3727 A, shifted around 8200 A for 
z~1.2 sources). We use spectra provided by the Keck DEEP2 
survey, obtained with the high spectral resolution DEIMOS 
spectrograph (Faber eFal. 20021) and taken along slits of typ- 
ical size 5" x 1". The data was reduced using a dedicated 
pipeline (Newman et al. 2012, Coope r et al. 2012| ). A velocity 
resolution of about 50 km s" 1 (R = A/AA « 5000) enables 
to resolve the [Oil] line (IDavis et al. 2003b . DEEP2 spectra 
are not flux calibrated and we use CFHT I band magnitudes 
(CFH1 2K camera) to c alibrate the galaxy-integrated "ID" spec- 
trum (ICoil et al. 2004b . 

Even if the [Oil] forbidden line luminosity is not directly 
coupled to the ionizing luminosity, it is possible to establish it 
empirically as a quantitative SFR tracer for the whole galaxy, 
and to express the SFR as 



SFR 

M Q .yr- 



= (5.3 ± 1.5) 10- 



-[OII] 



where L[ ii] is the [Oil] line luminosity calculated from the 
galaxy-integrated "ID" spectrum (Kennicutt 1998b]). 

The values obtained for the SFR have to be corrected for dust 
extinction. The Galactic extinction is taken into account in the 
C FHT magnitu des, as they already include the dust corrections 
of |Schlegel, Finkbeiner & Davis 19 98 , but the dust present in the 
galaxies themselves is an even larger cause of extinction. The 
above relationship between SFR and L[ ii] was directly obtained 
from a similar Ha calibrated relation, so has onl y to be corrected 
with the Ha extinction, A Htt (Kennicutt 1998b). The dust distri- 
bution was derived from the SED of the galaxies, obtained from 
the B, R and I band CFHT magnitudes with the kcorrect soft- 
ware developed by Blanton & Roweis 2007 (version v4_2). The 
software provides SED curves with and without extinction, from 
which we can easily deduce An a . 
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Fig. 2. [Oil] line and CO luminosities (respectively left and right panels for each galaxy) in position- velocity planes corresponding 
to the DEEP2 slits. Smoothed ensembles of clumps are separated by eye along the vertical axis, as shown with the white horizontal 
lines. The [Oil] diagrams were normalized in order to have fluxes proportional to the SFR, but all galaxies share the same arbitrary 
color scales. One arcsecond corresponds approximately to 8.5 kpc (Table [T] in the Appendix). 



The obtained extinctions are coherent with predicted 
extinction ratios dO'Donnell 1994b and observations (e.g. 
Forster Schre iber et al. 201TT) . but the inferred SFR (Table 
[T] in A ppendix) differ up to 3cr from the values ob- 
tained by ITacconi et al. 20131 These latter values derived from 
UV and Spitzer 24 jim luminosities, with the methods of 
Wu yts et al. 2011a give the same orders of magnitude, but 
have a larger scatter for the four galaxies considered here. The 
emp irical calibration t hrough B-band luminosity presented by 
Moustak as et al. 2 006 gives different values as well, the three 
methods being subject to comparable uncertainties and none of 
them being clearly better than th e others. The SFR obt ained by 
ITacconi et al. 20131 and through Moustakas et al. 2006 B-band 
calibration are given for information in Table [T] (Appendix), but 
would both increase the scatter in the resulting sub-galactic KS 
diagram without significantly changing our conclusions. 



From the CO and [Oil] lines, we estimate the gas mass and 
the SFR contained in areas of 1" in diameter (corresponding to 
the width of the slice and approximately to the size of the ensem- 
bles of clumps), and obtain the corresponding averaged surface 
densities, S gas and Ssfr- The resulting KS diagram is displayed 
in Figure [3] The data points scatter around the line of slope N=l 
corresponding to a depletion time (td ep i = M gas /SFR) of 1.0 Gyr, 
but this latter quantity is locally very different from one data 
point to the other. As shown in Figure|4] the mean depletion time 
is 1.4 Gyr, with a standard deviation of 0. 9 Gyr, as compare d for 
whole galaxies to 2.0+0.8 G yr at z=() in|Bigiel et al. 2008| 0.5- 
1.5 Gyr f rom z~2 to z~0 inlGenzel et al. 20101 and -0.7 Gvr at 
z=l-3 in ITacconi et al. 20T31 Our sample is too incomplete to 
compute a best fit slope, but the method developed here could be 
applied to more high redshift galaxies. 



3. A resolved Kennicutt-Schmidt law 

HST images of Figure [T] reveal kpc-sized clumps within diffuse 
regions, but these clumpy features (~ 0.1") are smoothed out at 
DEEP2 and IRAM resolutions, whose ranges are respectively 
0.6-1.0" and 0.5-1.6". Spectroscopy, however, helps to separate 
different components, thanks to their kinematics. DEEP2 spectra 
correspond to positio n- velocity diagram s (PV diagrams) along 
the galaxy major axis jDavis et al. 20071) . Figure [2] compares the 
[Oil] position-velocity diagrams with the corresponding slices 
in CO(3-2), both with the same 1" width. Smoothed 1 "-sized en- 
sembles of clumps are separated by eye along the velocity axis 
of the PV diagrams. The characteristics of these ensembles are 
given in Table 2 (Appendix), and we aim at describing the star 
formation efficiency within them. We tried to compensate the 
substructure separation by eye by taking into account all identi- 
fiable ensembles of clumps. 



4. Discussion and conclusion 

4.1. Advantages of the method 

We have shown how the various ensembles of clumps can be 
separated by their kinematics in PV diagrams, even though the 
angular resolution of molecular gas and SFR data was not able 
to separate them in integrated intensity. As such, the KS diagram 
can be obtained within regions of the resolution size (~1"). 

Previous resolved KS work at high redshift was only ob- 
tained using serend ipitous amplification by gravitational lenses. 
IDecarli et al. 2012I carried out one of the first spatially resolved 
studies in high redshift galaxies, using [Nil], FIR and CO obser- 
vations for two gravitationally lensed z~3.9 galaxies. They ob- 
tain a steep relation of slope N = 1.4 ±0.2 between the dust con- 
tinuum and the molecular gas surface brightness. Strong lenses 
are rare, and the determination of the clumps physical parame- 
ters depend on the lensing model. Our method is probably more 
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Fig. 3. Kennicutt-Schmidt diagram for 1" ensembles of clumps 
of EGS13003805 (green), EGS13004291 (red), EGS12007881 
(blue) and EGS 13019128 (purple), also displaying data points 
for these four galaxies as a whole (grey). The dotted diagonal 
lines correspond to constant gas depletion times of 0.1, 1, and 
10 Gyr from top to bottom, and the solid black line to the best 
linear fit with an N=l slope for the sub-galactic features. The 
corresponding depletion time is td ep i=1.0 Gyr. The error bars of 
0.3 dex yield a reduced^ 2 close to L The grey data points are 
taken from Gen zel et al. 2010l and lTacconi et al. 20131 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



t dep i [Gyr] 

Fig. 4. Distribution of depletion time scales for the 16 regions 
identified in the four galaxies studied in this letter. The distribu- 
tion has a mean of td ep i = 1 .4 Gyr and a standard deviation of 
0.9 Gyr. 

appropriate for a systematic study of the star formation at high 
redshift, until higher resolution instruments resolve the clumps. 

In the absence of high-resolution molecular gas data, 
Swinban k et al. 20121 report adaptive optics Ha observations of 
eleven kpc-scale star forming regions identified in z = 0.84-2.23 
galaxies and measure the velocity dispersion of the ionized gas <x 
and the star formation surface density 2 S fr. By assuming that cr 
also corresponds to the dispersion of the cold clumps and that the 
clumps are marginally stable with a Toomre parameter Q - h 
they claim a correlation between the gas surface density S gas and 
Ssfr- But the method is highly indirect, relies on many assump- 
tions, and underestimates beam smearing effects in the determi- 
nation of cr. 

4.2. Biases and uncertainties 

The four galaxies studied here are particularly luminous, and 
were selected from a sample dTacconi et al. 20101) with stel- 
lar masses and star formation rates respectively higher than 
3. 10 10 M© and 40 M yr" 1 . A high luminosity was also needed 



to visualize and isolate the clumps in the PV diagrams, so this 
method is intrinsically biased towards massive galaxies. 

The main uncertainties to estimate gas and SFR surface den- 
sities come from the values used for a = Mh 2 /L^ , for the CO(3- 
2)/CO(l-0) luminosity ratio, and for the extinction Ah<*. These 
quantities could vary significantly from one galaxy to the other, 
and would cumulatively infer multiplicative factors as high as 2 
or 3 in our final values for the mass of gas and of the SFR. 

4.3. Conclusion 

Our results, as well as most other observations 
(IGenzel et al. 20101 ITCcconi et al. 20T0l ITacconi et al. 20131 
IDecarli et al. 20121 ISwinbank et al. 2012b . indicate that the star 
formation scaling law between SFR and gas surface densities 
is not significantly different at high redshift than in the local 
Universe. Our limited sample of ~8 kpc-scale ensembles of 
clumps of distant galaxies is compatible with a power law 
S S fr °f exponent N=l and a depletion time of 1.0 Gyr, 

which is of the same order of magnitude as measurements 
at lower redshift. This adds to the growing evidence that the 
star formation processes ten billions years ago were similar to 
the ones that are observed in the local Universe. The method 
developed here could be applied to a more significant sample of 
high redshift galaxies, to obtain more statistically robust results. 
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Table 1. Some global properties and observational details for the four galaxies studied in this letter: (1) the redshift, as determined 
in the AEGIS catalogues from the [Oil] line using DEEP2 spectra ; (2) the angular distance Da and the luminosity distance D L ; 
(3) the gas mass derived from the IRAM Plateau de Bure observations ; (4) the SFR from the [Oil] line luminosity in the galaxy- 
integrated "ID" spectrum, obtained according t o section |2.2| ; (5) the tot al star mass obtained by SED fittings from CFHT B, R, 
and I bands photometric data using kcorrect JBlanton & Roweis"20 07). corrected for a Universe with Ho = 70 km. s -1 . Mpc -1 
; (6) the gas fraction f gas = M gas /(M gas + M star ) ; (7) the optical half light radius R\/2 as indicated in ITacconi et al. 20T31 ; (8) 
the surface densities associated to the gas mass, the star formation rate and the star mass, calculated within R\/2, for example 
as S gas = 0.5 M g3LS /nRy 2 ; (9) Ha, visible and [Oil] extinctions as derived from the kcorrect SED reconstruction ; (10) the 
extinction-corrected SFR obtained from 24 fim and UV continuum by ITacconi et al. 20"T3l ; and (11) the SFR ij/ extrapolated from 
the empirical calibration of Moustaka s et al. 20 06 (section 4 .3). A v , A[ n], SFR 2 4 yUm+ uv and ifs are given for information only. The 
comparison between the SFR obtained acco rding to section 2.2 SFR^m+uv and if/ gives an idea of the uncertainties in the SFR. 
Using SFR24 yUm +uv from ITacconi et al . 20 1 3 1 would give a constant depletion ti me linear fit with tdepi = 1.1 Gyr and values for td ep i 
between 0.6 and 10.3 Gyr (mean 2.1 Gyr), whereas using if/ calibrated through Moustakas et al. 2006l would give a fit with t^epi = 1.2 
Gyr and tdepi between 1.6 and 4.2 Gyr (mean 1.8 Gyr). 
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Table 2. Some properties of the clumps of EGS 13004291, EGS 13003805, EGS 12007881 and EGS 13019128 obtained from measurements in 
the position- velocity diagrams: (1) the gas mass obtained as a fraction of the total mass of the galaxy from the CO position- velocity diagram: (2) 
the SFR obtained similarly from the [Oil] position-velocity diagram, normalized with the total SFR determined according to section [2^2]; (3) the 
derived gas mass and SFR surface densities averaged over the same area of 1" in diameter ; (4) the depletion time t dep i = M gas /SFR ; and (5) the 
FWHM obtained from the IRAM position- velocity diagram with a gaussian fit in the direction of the slice (for information only). The FWHM 
roughly corresponds to the IRAM beam size, so we preferred to use the same 1" size for all clumps in the different calculations. Using the FWHM 
instead of a constant 1" diameter to calculate the surface densities gives a much larger scatter, but actually does not change the constant depletion 
time fit and the values of t dep i - The clumps are numbered from bottom to top according to the horizontal separation lines of Figure [2] 



Clump 


M gas [10 10 Mo] (1) 


SFR [M .yr _1 ] (2) 


/ ^gas 


(3) 


, / ^SFR 


(3) 


tdepi [Gyr] < 4) 


FWHM ["] (5) 


1Ogl0 lM G .pc-J 


1Og 10 Lr 11-2 

VMo.yr-^kpc l j 


EGS 13003805-1 


11.85 


51.17 


^ 01 Q 
j.yj iy 


— U. jM-D 


2.32 


0.93 


EGS 13003805-11 


6.92 


46.52 


2.785 


-0.387 


1.49 


0.78 


EGS 13003805-III 


3.18 


33.45 


2.447 


-0.531 


0.95 


0.63 


EGS 13004291-1 


2.37 


15.15 


2.324 


-0.870 


1.56 


0.78 


EGS 13004291-11 


5.80 


36.10 


2.713 


-0.494 


1.61 


0.64 


EGS 13004291 -III 


7.75 


34.88 


2.838 


-0.508 


2.22 


0.74 


EGS 13004291 -IV 


4.96 


31.89 


2.644 


-0.547 


1.55 


0.79 


EGS 13004291-V 


6.28 


45.10 


2.761 


-0.397 


1.44 


0.61 


EGS 13004291-VI 


3.99 


31.26 


2.550 


-0.556 


1.28 


1.1 


EGS 12007881-1 


5.33 


34.94 


2.680 


-0.503 


1.52 


1.35 


EGS 12007881-11 


7.26 


57.64 


2.814 


-0.286 


1.26 


1.36 


EGS 12007881-m 


1.41 


70.27 


2.10 


-0.200 


0.20 


1.57 


EGS 13019128-1 


1.70 


5.43 


2.165 


-1.330 


3.12 


0.83 


EGS 13019128-11 


0.58 


77.22 


1.697 


-0.177 


0.07 


1.14 


EGS 13019128-m 


7.08 


71.22 


2.785 


-0.212 


0.99 


1.67 


EGS 13019128-IV 


5.25 


64.10 


2.655 


-0.258 


0.82 


1.60 


EGS 13019128-V 


0.78 


43.61 


1.828 


-0.425 


0.18 


1.62 


EGS 13019128- VI 


2.60 


7.94 


2.351 


-1.165 


3.28 


1.86 



